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Abstract 

Integrable quantum mechanical systems with magnetic fields are 
constructed in two-dimensional Euclidean space. The integral of mo- 
tion is assumed to be a first or second order Hermitian operator. Con- 
trary to the case of purely scalar potentials, quadratic integrability 
does not imply the separation of variables in the Schrodinger equation. 
Moreover, quantum and classical integrable systems do not necessarily 
coincide: the Hamiltonian can depend on the Planck constant h in a 
nontrivial manner. 
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1 Introduction 



The purpose of this article is to study the integrability properties of a quan- 
tum particle moving in an external magnetic field. More specifically, we will 
consider the Schrodinger equation in a two-dimensional Euclidean space with 
the Hamiltonian 

H = -^(d 2 x +d?)-^[A(x, y)d x +d x A(x, y)+B{x, y)d y +d y B(x, y)]+V(x, y). 

(1.1) 

The vector and scalar potentials (A, B) and V are to be determined from 
the requirement that the system should be integrable, i.e. a well-defined 
quantum mechanical operator X should exist, that commutes with the Hamil- 
tonian, i.e. 

[H,X]=0. (1.2) 

In this particular study, we shall restrict to the case when X is a first or 
second order polynomial in the momenta. We shall be particularly interested 
in the case of superintegrable systems, when two independent operators, X\ 
and X 2 , commuting with the Hamiltonian exist. In general, X% and X 2 do 
not commute with each other, but together generate an algebra of operators, 
commuting with H. 

In classical mechanics, integrable systems are of interest, because they 
have regular trajectories. Indeed, their motion is restricted to a torus in 
phase space. Superintegrable systems are even more regular. Trajectories 
are completely determinded by the values of the 2n — 1 integrals of motion. 
In particular, all bounded trajectories are periodic, as in the case of the 
harmonic oscillator, or Kepler problem. 

In quantum mechanics, integrability, i.e. the existence of n integrals of 
motion, provides a complete set of quantum numbers, characterizing the sys- 
tem. Moreover, it simplifies the calculation of energy levels and wave func- 
tions. Superintegrability, in all cases studied so for, entails exact solvability. 
This means that energy levels in superintegrable systems can be calculated 
algebraicly, i.e. they satisfy algebraic rather than transcendental equations. 

Previous searches for integrable and superintegrable systems in quantum 
mechanics concentrated on scalar potentials only (Tj-IIIlj. It was established 
that for scalar potentials the existence of first and second order integrals 
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of motion implies the separation variables in the Schrodinger equation, and 
also in the Hamilton- Jacobi equation in classical mechanics. Moreover, for 
scalar potentials and second order integrals of motion, classical and quantum 
integrable systems coincide (i.e. classical and quantum potentials are the 
same). 

Surprisingly, when third order integrals are considered, a new phenomenon 
occurs: integrable and superintegrable quantum systems that have no clas- 
sical counterpart [H]-[Ilj. Indeed, in the classical limit ft — > the potential 
vanishes, V(x, y) — > and we obtain free motion. 

Previous studies of integrability in magnetic fields were conducted in the 
framework of classical mechanics [Hp, [IE|. It was established that the exis- 
tence of second order integrals of motion in the presence of magnetic fields no 
longer implies the separation of variables. However, the integrals of motion 
were still classified into equivalence classes under the action of the Euclidean 
group and the highest order terms have the same form as in the case of a 
purely scalar potential. 

In this paper we restrict ourselves to the two-dimensional Euclidean space 
E(2), the Hamiltonian Ijl.lj) and to first, or second order integrals. In Sec- 
tion 2 we formulate the problem of finding the integrals of motion, first in 
the classical, then in the quantum case. We show that the determining equa- 
tions in the two cases are the same for first order integrals of motion, not 
however for second order ones. Section 3 is devoted to first order integrals 
of motion. They are shown to exist if and only if the magnetic field and an 
effective scalar potential are invariant under either translations, or rotations. 
We also show that superintegrability with two (or more) first order integrals 
occurs only for a constant magnetic field and effective potential. In Section 
4 we consider a specific class of second order operators which we call "carte- 
sian integrals". In the absence of a magnetic field they lead to separation of 
variables in cartesian coordinates. We also show that superintegrability with 
one cartesian integral and a second integral of any (quadratic) type occurs 
only for a constant magnetic field. In the cartesian case there is no difference 
between classical and quantum integrability. Polar integrability and super- 
integrability are investigated in Section 5. All cases of integrability with one 
"polar" integral of motion are identified. The quantum case differs from the 
classical one and the magnetic field can depend on the Planck constant % in 
a nontrivial manner. In Section 6 we show that a polar integral can exist 
simultaneously with any other independent second order integral only if the 
magnetic field is constant. The final Section 7 is devoted to conclusions and 
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open problems. 



2 Formulation of the problem 
2.1 Classical mechanics 

Since we will be comparing results in quantum and classical mechanis, let 
us briefly recapitulate some results obtained earlier (THJ , [IE)- The classical 
counterpart of the Hamiltonian (II. lj) is 

H = \{vl + P 2 y ) + A(x, y)p x + B{x, y)p y + V(x, y), (2.1) 

where p x and p y are the momenta conjugate to x and y, respectively. The 
classical equations of motion in the Hamiltonian form are 



*=% = p*+a, y = ^ = p y + B 

PX Q x A X P X B X Py, Py Vy AyP X ByPy 

(2.2) 

The equation of motion (12 .2J) can be rewritten in the Newton form as 

■± = -W x + Qy, , , 

y = -Wy - fix, { ' 

W = V-l(A* + Bi), {2A) 

fi = Ay — B X . 

The equations of motion (12 .'M are invariant under a gauge transformation 
of the potentials 

V{x, y) -> V(x, y) = V + (A V 0) + |(V0) 2 , (2 g) 

where we have put A = (A, B) and = (f)(x, y) is an arbitrary smooth 
function. Thus, the quantities that are of actual physical importance are the 
magnetic field Q and the effective potential W. 

A classical first integral of motion is postulated to have the form 

C = fi(x,y)x + f 2 (x,y)y + m(x,y). (2.6) 
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The determining equations for the funtions fx, f 2 and m are obtained 
from the requirement 




(2.7) 



when eq. ()2.2|) are satisfied, i.e. C is a constant on the solutions of the 
equations of motion and Poisson commutes with the Hamiltonian. 
Similarly, a classical second order integral of motion has the form 

C = gi(x, y)x 2 +g 2 (x, y)y 2 +g 3 (x, y)xy+k 1 (x, y)x+k 2 (x, y)y+m(x, y). (2.8) 

The determining equations for the functions ki and m are again ob- 
tained from the condition (12. 7j) . 

The equations for the coefficients of the first and second order classical 
integrals of motion were derived and partially solved elsewhere [11], [Tl]. We 
shall give them again below as classical limits of the corresponding equations 
in the quantum case. To facilitate a comparison, we must rewrite the classical 
integrals in terms of momenta, rather than velocities, i.e. substitute x = 
p x + A, y = p y + B. 

2.2 Quantum mechanics 

In quantum mechanics an integral of motion will be a Hermitian operator X 
that commutes with the Hamiltonian H. 

Let us first consider a first order integral in the momenta: 



The classical limit of the operator (|2.9J) is the integral (|2.6j) ; fx, f 2 and m 
are functions of x and y. 

The commutator [X, H] with H as in eq. (|l.ll) will contain second, first 
and zero order terms in the derivatives. Setting the coefficients of all of them 
equal to zero, we obtain the following set of determining equations 



X = -y (fid x + d x fx + hd y + d v f 2 ) + fxA + f 2 B + m. 



(2.9) 



fir 



0, f2,y — 0, fx,y + f2,x — 0, 



(2.10) 



-f 2 Q + m x = 0, fxQ + m y = 0, fxW x + f 2 W y = 0. 



(2.11) 
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We see that the Planck constant fi does not figure in eq. (I2.1()|) and (|2.11j) . 
Hence these equations must coincide with their classical limit, and indeed, 
they do [TH]. In particular, eq. (|2.10|) implies 

f = ay + (3, f 1 = -ax + 1 , (2.12) 

where a, (3 and 7 are real constants. Hence the leading terms (independent 
of A, B and m) of the operator X of eq. (|2.9|) lie in the Lie algebra e(2) of 
the Euclidean group E(2), generated by 

Pi = -i?id x , P 2 = -ihdy, L 3 = -ih(yd x - xd y ). (2.13) 
Thus, we have 

X = aL 3 + (]P 1 + jP 2 + a(yA - xB) + @A + jB + m. (2.14) 

We shall write the second order operator corresponding to the integral (|2.8j) , 
after symmetrization, as 



X = -\h 2 {2 9l d 2 x + 2g 2 d 2 + 2g 3 d x d v 
+(2#i,z + g3, y )d x + (2g 2 , y + g 3 ,x)d y + gx >xx + g 2<yy + g 3<xy } 
-%{{4giA + 2g 3 B + 2k x )d x + {Ag 2 B + 2g 3 A + 2k 2 )d y (2.15) 
+2gtA x + 2g hx A + 2g 2 B y + 2g 2 ^ y B + g 3 A y + g 3 B x + Ag 3 , y + Bg 3jX 
+h, x + k 2 , y } + gi A 2 + g 2 B 2 + g 3 AB + k ± A + k 2 B + m 

The commutativity condition [H, X] = implies the following set of 
determining relations: 

9i,x = 0, g 2 , y = 0, g hy + g 3jX = 0, g 2iX + g 3)V = 0, (2.16) 



hx ~ 93^ = 0, k 2y + g 3 £l = 0, 

2n(g 1 - g 2 ) + ki y + k 2x = 0, , . 

2giW x + g 3 W y + k 2 n -m x = 0, {Z - U > 
+ g 3 W x - hQ -m y = 0, 

h 2 

hW x + k 2 W y + —(g 2x n y - g ly n x ) = 0. (2.18) 

Eq. (12.161) and (I2.17jl are the same as the classical ones [11]. Eq. (|2.18j) is 
however different. It involves the Planck constant and reduces to the classical 
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case only in the limit fi — > 0. Thus, in the presence of a nonconstant magnetic 
field £l(x,y), classical and quantum integrability differ! 

Eq. ()2.16|) can be solved as in the classical case [TH] and they imply 

91 = ®y 2 -Py + 8, 

g 2 = ax 2 + 72 + C, (2.19) 
g 3 = -2axy + (3x - jy + f , 

where the Greek letters represent real constants. Substituting (|2.19jl into (|2.15j) 
we obtain the operator X in the form 

X = a[(L 3 + yA- xB) 2 + Ti 2 } - \(3[{L Z + yA- xB)(P 1 + A) 
+(P 1 + a) (L3 +yA- xB)] - ± 7 [(L 3 + yA - xB)(P 2 + B) + (P 2 + B) 
(L 3 + - xS)] + 5(Pi + A) 2 + C(P 2 + P) 2 + ^(Pi + A)(P 2 + 5) 



f (2hd x + k l)X + 2k 2 d y + k 2 , y ) + kiA + k 2 B 



m. 



(2.20) 

Thus, the leading part of eq. ()2.20jl lies in the envelopping algebra of e(2). 
For A = B = this coincides with the case of a scalar potential [I], 0, [T7] . 
As in the scalar case we can simplify eq. (|2.2()l) by Euclidean transformations 
and linear combinations with the Hamiltonian. The operator X is trans- 
formed into a similar operator, with new values of the constants a, 
Four classes of such operators exist, represented by 



X c = (Pi + A) 2 -f(2k 1 d x + k hx 
+2k 2 d y + k 2 , y ) + hA + k 2 B + m, 



(2.21) 



X R = (L 3 +yA- xB) 2 + h 2 - f(2hd x 

+h, x + 2k 2 d y + k 2>y ) + hA + k 2 B + m, [ ' > 

X P = -I{(L 3 + yA - xB){P x + A) + (P, + A)(L 3 + yA - xB)} 

-f{2hd x + k ljX + 2k 2 d y + k 2 , y ) + hA + k 2 B + m, { 6) 



X E = (L 3 + yA - xB) 2 + h 2 + a[(Px + A) 2 - (P 2 + P) 2 ] 
-f (2&A + k liX + 2k 2 d y + k 2)V ) + k x A + k 2 B + m, a> 0. 



(2.24) 



In the case of a purely scalar potential the existence of a commuting op- 
erator of the type Xq, Xp, Xp or Xp implies that the Schrodinger equation 
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will allow separation of variables in cartesian, polar, parabolic or elliptic co- 
ordinates, respectively. In the last case a is related to the interfocal distance 
for the elliptic coordinates. 

Substituting eq. (|2.19ft into eq. 1)2. 18j) we obtain 



Thus, for the special case a = (3 = 7 = classical and quantum integra- 
bility will coincide. 



3 First order integrability and superintegrabil- 



A first order integral (|2.9J) in quantum mechanics will exist if the overde- 
termined system ()2.1()l) and (12.11)) has a solution. The general solution of 
eq. (I2.1()j) is given by eq. ()2.12l) . Substituting into ()2.11l) we obtain: 

(ax - 7)^ + m x = 0, (ay + + m y = 0, (ay + [3)W X + (-ax + j)W y = 0. 

(3.1) 

We are only interested in cases with a magnetic field present, i.e. Q ^ 0. 
With no loss of generality, we need only distinguish two cases: 
1. a = 1, P = j = 0. 
We obtain 



We see that in this case the magnetic field Q and the effective scalar poten- 
tial must be spherically symmetric. The potentials in the Hamiltonian (jl.ll) 
can by a gauge transformation be taken into 




(2.25) 



ity 



w = w( P ), n = Q( P ), P = v^ 2 + y 2 - 



(3.2) 




(3.3) 



The integral of motion is 




(3.4) 
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2. a = /3 = 0, 7 = 1. 

The magnetic field and effective potential are translationally invariant 



Q = Q{x), W = W{x), (3.5) 

and we can take 

A = Qy, B = 0, V = W + \y 2 tt\ 
X = P 2 + jQ{x)dx. 

The system Ijl.lj) will be first order superintegrable if at least two first 
order integrals (|2.14jl exist. This is only possible if the magnetic field and 
effective potential are constant: 

ft = n , W = W . (3.7) 

In this case actually three operators commuting with the Hamiltonian 
exist and we have 

H = ±(ihd x -ny) 2 -¥-dl (3.8) 

X 1 = P u X 2 = P 2 + xn , X 3 = L 3 - i(x 2 - y 2 )n , (3.9) 
where we have chosen the gauge to be such that 

A = Q y, B = 0, V=hl 2 y\ (3-10) 



The classical equations of motion (|2.2j) are easily solved. The trajectories 
are circles (and are hence all closed). The Schrodinger equation allows the 
separation of variables in cartesian coordinates. The solution is 

y(x,y) = e l ^f(y), (3.11) 
where f(y) satisfies the harmonic oscillator equation 



/"-{(^) 2 -f}/ = 0. (3.12) 
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The integrals of motion (J3.9I) . together with the constant f2, satisfy the 
commutation relations of a central extension of the Euclidean Lie algebra: 

[X U X 2 ] = -iMh, [X 3 ,Xx] = -ihX 2 , [X 3 ,X 2 ] = ihXx. (3.13) 

Only three of the integrals X%, X 2 , X 3 and H can be independent and 
indeed they satisfy 

X\ + X\ + 2fi ^3 -2/7 = 0. (3.14) 
In polar coordinates the Schrodinger equation 

{|(^cos(0)<9 r - S^fy - fir sin(0)) 2 

-i( S in(0)9 r + ^%) 2 }* = ^ 1 j 

R-separates [12], [TH], rather than separates, and we have 

qr(r, 0) = e-i° r2sin2 ^J m (A;r)e fm ^, A; 2 = IE + mfi, (3.16) 
where J m (kr) is a Bessel function. 

4 Cartesian integrability and superintegrabil- 
ity 

4.1 Integrability 

In order to find integrable systems with a second order operator commuting 
with the Hamiltonian, we must solve the system (I2.16|) to (I2.18J) . To do this, 
we first transform X to its canonical form, i.e. one of ()2.2ip to (j2.24ji . We 
start with the simplest case, namely X c of (|2.21|) . We call this the "cartesian" 
case, because for a purely scalar potential it corresponds to separation of 
variables in cartesian coordinates. It corresponds toa = /? = 7 = C = 
£ = and 5 = 1 in eq. (|2.20jl . Eq. (|2.25j) implies that the determining 
equations (|2.16|) , (j2. 17j) and (|2.18j) are the same in the classical and quantum 
cases (the h 2 term in eq. 1)2. 18j) vanishes). For purely scalar potentials fi = 
k\ = k 2 = we reobtain the known result W = W (y) +m(x) [lj. From now 
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on we assume Q ^ 0. For completeness, we reproduce the result obtained 
earlier jTH] in the classical case, since it is valid in the quantum well: 



W 




m = -§Q? 3 + 2f - 3gf) + b(fg - / 2 ) + f (g 2 - / 2 ) + c(g - 2f) + eg - kf. 



Here a, 6, c, d, e and are constants and the functions / = f(x) and 
g = g(y) satisfy 



Two exceptional cases occur when we have f x = or g y = 0. These 
however imply Q = fl(x), W = W(x) or Q = £l(y), W = W(y) respectively. 
Then a first order invariant exists and the second order one is simply its 
square. The general solution of eq. (|4.2jl are elliptic functions. 

4.2 Cartesian superintegrability 

We shall now assume that Q and W are such that one cartesian integral X\ 
exists, i.e. they satisfy eq. (|4.1|) . We require that a second integral X 2 of the 
type (I2.2()|) should exist, in addition to the considered cartesian one. We can 
simplify the integral X 2 by translation and by linear combinations with Xi 
and H. Rotations can not be used, since they would change the form of the 
operator X\ and of the Hamiltonian. Two cases must be considered, a^O 
and a = 0. 

Case 1: a^O 

We set a = 1, by a translation we transform (/?, 7) — > (0,0), by linear 
combinations we set (5, C) "~ * (0,0). We are left with an operator X 2 in the 
form (|2.20|) with a = l,/? = 7 = 5 = £ = 0. The constant £ and functions 
fci, /C2 and m must be determined from the system (I2.17L (I2.18jl . Let us 
consider the case when Q and W are as in eq. (I4.1|h The first two equations 
imply 



(4.1) 



XX 



af 2 + bf + c, = -ag 2 + dg + e, 

/a 7^ 0, ^ ^ 0. 



(4.2) 



fci = -2y(xf x - /) - x 2 ?/^ + £/ a . + £x£ w + Ci(y) 
^2 = xy 2 fxx + 2x(yg y - g) - iyf xx - ig y + C 2 (x). 



(4.3) 
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We substitute k\, k 2 , £1 and W into the remaining four equations and in- 
vestigate their compatibility. After somewhat lengthy calculations we obtain 
a simple result: the equations are compatible for Q ^ if and only if Q and 
W are constant. We arrive at the case (3.7), already investigated in section 
3. 

Case 2: a = 

In order to obtain an independent second order integral we must have 
P 2 + 7 2 7^ and we can normalize j3 2 + j 2 = 1 and put 5 = £ = £ = (by 
linear combinations with H and X\). The set of equations 1)2.16)1 to 1)2.18)1 is 
then again compatible only for Q and W constant. 

The conclusion of this section is that for Q ^ cartesian superintegra- 
bility with two second order integrals exists only in a trivial sense. Thus 
Q = Qoi W = Wo and all second order integrals are reducible: they are 
polynomials in the three first order ones. 

5 Polar integrability 

We now request that one second order integral should exist and that it be 
of the form 1)2.22(1 . We shall call this operator X R a polar type integral. Let 
us transform the determining equations 1)2. 17)1 . 1)2.18)1 to polar coordinates 
x = rcos(0), y = rsin(0). The resulting equations are 

P r = 0, P + Q^ = 0, (5.1) 

2r 3 Vl -Pf- rQ r + Q = 0, (5.2) 

m - 2r 2 W^ + rPVL = 0, m r - QVt = 0, (5.3) 
h 2 1 

—% + PW r + -QW+ = 0, (5.4) 

where we have put P — k\ cos(0) + k 2 sin(^) and Q = — k\ sin(0) + k 2 cos(</>). 

We see that eq. 1)5.4)1 contains a term proportional to fi 2 . It follows that in 
this case quantum integrable systems will differ from classical ones, at least 
if we have ^ 0. In the classical limit ft — > the quantum systems will 
reduce to classical ones, or to free motion. This is a new phenomenon. In 
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the absence of magnetic fields, classical and quantum systems with second 
order integrals of motion coincide. 
Eq. 1)5.1)) imply 



P = -/'(</>), Q = f((l>) + R(r), (5.5) 

with /(</>) and R{r) to be determined. We shall use primes and dots to denote 
derivatives with respect to <ft and r, respectively. We again assume Q ^ 0. 
Indeed, for Q = we obtain the known case of a scalar potential, separable 
in polar coordinates: W = Wo(r) + Wx r 2 . 
We solve eq. 1)5. 2 j) for the magnetic field 

V = -l^(f" + f + R-rR). (5.6) 
From eq. ()5.3j) we obtain a compatibility condition (m r ^ = m^ r ), namely 



2r 2 ^+W + ^(r + /+i?-ri?) + ^ + ^(/ +J R)(r + /') =0. (5.7) 

Using eq. 1)5 .4jl to eliminate from eq. 1)5.7)1 . we obtain, for 

/ + ^0,/'/0 (5.8) 
an equation for W(r, <fi) that we can solve 

(3(f+R)-rR) w/ _ h 2 R(f'"+f ) , 3(/+fl) / r;/ , / , d , d\ , 



(5.9) 



We obtain 



VV ~ 8r 2 /' 32r 4 8^' ' U 8^ U + °' "V 32/'r 4 f5 10) 

-§£ + (^ + ^)ru + ^u- 1 [;(r 3 ii + 3r 2 u) + Fr 3 u + W , 

where F = F(0) and W (4>) are two new functions, introduced as integration 
"constants". We have also introduced the functions u{r) and S(r), satisfying 

S(r) = \R(r), u(r) = \s(r) (5.11) 
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Let us first consider the two special cases in eq. (15 .8j) . 
For /(0) + R(r) = eq. implies f2 = and we are not interested in 

this case. 

In the case /'(</>) = we have 

P = o, Q = Q(r), W = W(r), m = g, (g u) 

h = -Q(r)sin(0), A; 2 = Q(r)cos(0), 

and the classical and quantum cases coincide. Moreover, a first order integral 
exists. 

Let us return to the generic case ()5.10|) with conditions ()5.8|) satisfied. 
We substitute W of eq. (|H.10jl into eq. 1)5. 3j) and 1)5.4)1 to obtain 

m = ^ + i^ + ^- ^- + i^ + m o(^)' ( 5 - 13 ) 
m' = 2/'F, (5.14) 



... 6 r 3 u 2 J"\,l ,6 C 4F\ 3/'" 45 4tf' 

M + 2~(jr) J + ^ + ^ + jr) u + ^JJ u + ^FJJ + ^ + ^3j7 = °> 

(5-15) 

where 

= 9//"' _ f ( r_y _ 15/7" _ M£i 

32 32 \ /' / 32 4 ' 

B=f(f)'-if-^ ; (5.16) 

C = 6f + /(f)', 

The next task is to solve eq. 1)5.151) . Notice that this is not a partial 
differential equation. It involves four unknown functions u(r), f{4>), F(<j>) 
and Wo(4>), each depending on one variable only. Hence, we can consider 
this equation to be an ordinary differential equation for u(r), and then es- 
tablish the compatibility conditions on the other unknowns for which the 
dependence will cancel. The complete analysis is rather lengthy and involves 
the consideration of many special cases. We shall only present the main 
arguments and final results. 



Case 1: 



((y)'(^))' = 0, (5.17) 
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All subcases lead to the following solution (or special cases thereof): 

/ = C O + Cicos(0)+C 2 sin(0), , v 

F = K 1 , W = K 2 f + K 3 . ( °- La) 

Eq. (15. 15 J) then reduces to 



. 6 3. 3 15C 6K! 4K 2 

u + -u+ —u- —u- —- _ + __ = o. 5.19 

ty y>Z iy*0 Qf* fjf**J <yO 

The general solution of eq. (|5.19l) is 

C 2K X AK 2 be , 

8r 4 r l 3 r r 6 



and hence we have 



S = gfr - AK X + ar 3 - 6r - f 
= -C + 3ar 5 - 6r 3 + 3cr. 



(5.21) 



Finally, the magnetic field and effective potential in this case are 

Q = 6ar 2 - b, (5.22) 



W = -2ar(C\ cos(0) + C 2 sin(0)) + yr 4 - 3acr 2 - a 2 r 6 . (5.23) 

Since does not depend on the classical and quantum cases are the 
same. The corresponding classical integral of motion is 



Cr = (xy — yx) 2 + (—C2 — (3ar 5 — br 3 + 3cr) sin(0))i; 

+{Ct + (3ar 5 - br 3 + 3cr) cos((j)))y - ^ + 
_3^ri + 2C 1 ar 3 cos(0) - C x br cos(0) + 2C 2 ar 3 sin(0) ^ ^ 

-C 2 6r sin(0) + ^f- + ^ + ^ . 



Case 2: 



((y)'(]y))V0 (5.25) 



15 



A complete analysis [OE] shown that eq. (j5.15j) in this case is consistent 
only if we have 



and hence 



« = ^I + ^ + c (5.26) 



S = ~-b,R = a. (5.27) 



Moreover, the function /(</>) must satisfy 



£111 

(f + «) 2 (y )' + 24/'(/ + a) + 9f'(f + a) + 15/7" = 0. (5.28) 
The functions F(4>) and W ((j)) are given explicitly in terms of f(<f) as 

F=~^+ ^((f + «) 2 f - 2(/ + a)/" + iff) + ^ 

Wo = + ^ " T + C 2 - ( ' 

We integrate eq. (|5.28jl twice and put y = f + a to obtain the second 
order equation 

2, , N9 4A B 2 -A 2 , 

y" = —{y'f - 3?/ + — + — — 5.30 

2/ y y 3 

where A and I? are constants. 

This equation has a first integral K, in terms of which we have 

y\y'f = -y 6 + 2 V + (b 2 - A 2 ) y 2 + k. (5.31) 

This equation can be written as a quadrature that will express the in- 
dependent variable as a function of y in terms of elliptic integrals. The 
results are not very illuminating, so instead of presenting them, we restrict 
ourselves to some special cases. Let us first rewrite eq. (I5.31j) as 

y\y'f = -{y 2 - y\){y 2 - vlW - yl) = t{ v ), (5.32) 
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where the roots y\, y2 and y 3 are related to the constants A, B and K by the 
formulas 



K = ylylyl B 2 -A 2 = -{y\y\ + y\y\ + y\y\), 2A = y\ + y\ + y\. (5.33) 

If all the roots yi are real, the behavior of the polynomial T(y) as a 
function of y is shown on Fig. |l(a)[ 

If all roots are distinct (0 < y 3 < y 2 < y\ < oo), real periodic solutions are 
obtained for — y 3 < y < y 3 , y 2 < y < y\ and — y\ <y<y 2 . However, these 
are expressed in terms of elliptic functions and the period is not a multiple 
of 7r. Constant solutions of eq. ()5.32|) are obviously y = ±yk, k = 1, 2 or 3. 

Elementary dependent real finite periodic solutions are obtained when- 
ever the polynomial T(y) has multiple roots. The corresponding solutions 
are 

(1) 2/3 = 2/2 = 0, yxX) (See Fig, p^b^ 

y = yi sin(0 - O ) (5.34) 

(2) = t/ 3 <2/2<2/i (See Fig.|T^ 



y = ± ^yy 2 i + vl + (vl - vl) sin2 (^ - ^) ( 5 - 35 ) 

or in terms of A and B: 

y = ±^ A + B sin2((j) - cj) ) (5.36) 

(3) < y 3 < y 2 = Vl (See Fig. [1(d)]) 

In this case we give the solution y implicitly as 



-2^y\-y 2 arcsin(^) + yi(arcsin( ^±g^ 



-arcsin(J-5i T ) = ±2y^i( 



(5.37) 



The solution is real, finite and periodic for — y 3 < y < y 3 . 
For any solution y((f>) of eq. 1)5.31)) we obtain a magnetic field and effective 
potential in the form 

Q f" + f + a ^ 
= 2? ( } 
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rxr h " M 2 f" U'Y v {f + a) 2 J'" A , 3/"(/ + a) d 

W = — flH — — )—— — - — h4) — — \-Co 

8r 21 f + a {f + a) 2) 32r 4 1 /' 32r 4 2r 2 (/ + a) 2 2 



The functions P, Q and m figuring in the polar integral are 



(5.39) 



p = -f (0), Q = f(<j))+a, 

m = >ry -T -itSf + jctM 2 */ + " (/') 2 >- <5 ' 40) 

Let us sum up the results of this section. Three different cases of polar 
integr ability exist. They are given by eq. (ETH1 . (I5~2"2l to (I5~2"3I) and (IP£1) 
to (I5.4()jl . respectively. The last case provides an example where the quan- 
tum system and the classical one differ. Indeed, the Planck constant figures 
explicitly in the effective potential W and in the integral of motion. 



6 Polar superintegrability 

Let us assume that we have a Hamiltonian that is "polar integrable", 
i.e. allows an integral of motion of the form X R as in eq. 1)2.22)1 . The 
magnetic field Q and effective potential W must hence have one of the three 
forms established in Section 5. For the system to be superintegrable, it must 
allow at least one further integral, by assumption of the form 1)2.20)1 . We can 
simplify this second integral by linear combinations with Xr and with H and 
also by rotations, since they will not destroy the form of X R (nor H). Thus, 
in eq. 1)2.20)1 we take a = 0, ( = —5. Furthermore, we can assume (3 2 + r y 2 ^ 
0, since otherwise we would be in the case of cartesian superintegrability, 
already treated in Section 4. By a rotation and normalization, we can set 
P = 1, 7 = 0. It follows that the second integral X 2 is of the parabolic type, 
conjugate to X P of eq. 1)2.23)1 . 



The determining equations for X 2 , obtained from 1)2.17)1 and 1)2.18)1 are 

h, x ~{x + 0^ = 0, k 2>y + (x + 0^ = 0, (6.1) 
-2/3Qy + k hy + k %x = 0, (6.2) 



-2yW x + (x + g)Wy + k 2 tt -m x = 0, 
(x + £)W X - hQ - m y = 0, 



(6.3) 
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hW x + k 2 W y + —pQ x = 0. (6.4) 

For each of the three polar integrable systems we obtain the same result, 
namely: equations (|6.ip to ()6.4|) are compatible only if Q = f2 and W = Wo 
are constant. Then we have three first order integrals and the corresponding 
second order integrals are polynomial in the first order ones. 



7 Conclusions 

We have constructed all integrable quantum systems with a vector and scalar 
potential (as in eq. (II. lj) ) that possess either a first order integral, or a second 
order one of the cartesian, or polar type. 

It is interesting to compare such systems with a nonzero magnetic field 
fl with systems allowing a scalar potential only. 

1. The first difference is that for Q ^ quantum and classical integrable 
systems with second order integrals do not necessarily coincide. The Planck 
constant % can figure in a nontrivial way in the potentials and integrals of 
motion. 

2. The existence of a first order integral of motion implies a geometrical 
symmetry, both for Q ^ and Q = 0. Indeed, a first order integral exists if 
and only if we have either Q = Q(r), W = W(r), or Q = Q(y), W = W(y) 
(up to Euclidean transformations). The functions Q and W are arbitrary in 
both cases. 

3. The existence of a second order integral for Q = implies that the 
Schrodinger equation will allow separation of variables in cartesian, polar, 
parabolic, or elliptic coordinates. In each case the potential V(x, y) depends 
on two arbitrary functions of one variable. For Q ^ the coordinates no 
longer separate. The requirement that an irreducible second order integral 
should exist for ^ is much more restrictive than for Vt = 0. The quant- 
ities fl(x,y) and W(x, y) again depend on two functions of one variable, 
however these functions obey certain ordinary differential equations. They 
are hence determined completely, up to some arbitrary constants. For in- 
stance, in the cartesian case, they are elliptic functions, or degenerate cases 
of elliptic functions. 

4. For Q = four families of superintegrable systems in E(2) exist [l], 
each depending on three parameters. For O ^ we have shown that super- 
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integrability with first order integrals of the cartesian, or polar type, exists 
only for Q and W constant. 

Several related problems are presently under consideration. To complete 
the study of quadratic integrability in E{2) for fl ^ we must still consider 
parabolic and elliptic integrability. For Q = their is a close relation between 
superintegrability and exact solvability §251 • For Q ^ the requirement 
of superintegrability seems to be too restrictive. An important question is 
whether some of the integrable systems found in this article are actually 
exactly solvable. 
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(a) Three pairs of simple roots 



(b) One quadruple root, a pair of single 
ones 




(c) One double root, two pairs of single (d) One pair of double roots and one of 
ones simple ones 



Figure 1: Roots of the polynomial T(y) in eq. (|5.32l) 
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